source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/service_plot_gdp.R')
source('~/Desktop/service/figures/service_plot_gdp.R')
source('~/Desktop/service/figures/us_sector_plot.R')
source('~/Desktop/service/figures/us_sector_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
View(d)
source('~/Desktop/service/figures/us_intensity_plot.R')
View(d)
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
library(data.table)
library(here)
library(magrittr)
dir = here()
setwd(dir)
va = fread("real_va_1997.csv")
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/global_intensity_plot.R')
source('~/Desktop/service/figures/global_intensity_plot.R')
source('~/Desktop/service/figures/global_intensity_plot.R')
source('~/Desktop/service/figures/global_intensity_plot.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
source('~/Desktop/service/figures/sector_energy_plot.R')
source('~/Desktop/service/figures/sector_energy_plot.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
x = log(d$SL.SRV.EMPL.ZS)
y = log(d$EG.USE.PCAP.KG.OE)
r = lm(y ~ x)
r.square = summary(r)$r.squared
x = seq(log(1), log(100), length.out = 50)
x = data.frame(x)
regression = exp(predict(r, x,  interval = "prediction", level = 0.99))
regression = data.frame(x = exp(x), regression)
energy_total.r2 =  paste("R^2 == ", round(r.square,2))
source('~/Desktop/service/figures/absolute_dematerial.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
library(data.table)
library(ggplot2)
library(here)
library(gridExtra)
library(grid)
library(gtable)
library(reshape2)
text.size = 10
# global
###################################################
dir = here()
wd = gsub("figures", "data/world bank", dir)
setwd(wd)
d = fread("data_format.csv")
plot(d$NV.SRV.TETC.ZS, d$EG.EGY.PRIM.PP.KD, log="xy")
plot(d$NV.SRV.TETC.ZS, d$EG.EGY.PRIM.PP.KD, log="xy", xlim = c(0.01, 100))
plot(d$NV.SRV.TETC.ZS, d$EG.EGY.PRIM.PP.KD, log="xy", ylim = c(0.01, 100))
plot(d$NV.SRV.TETC.ZS, d$EG.EGY.PRIM.PP.KD, log="xy", ylim = c(0.1, 100))
plot(d$NV.SRV.TETC.ZS, d$EN.ATM.CO2E.PP.GD.KD, log="xy", ylim = c(0.1, 100))
source('~/Desktop/service/figures/relative_dematerial.R')
png("relative_dematerial.png", width = 7.5, height = 3.2,  units = 'in', res = 600)
grid.arrange( rbind(
cbind(gA, gB, size = "last"),
cbind(gC, gD, size = "last"),
size = "last")
)
dev.off()
# export
############################################################
setwd(dir)
gA = ggplotGrob(energy)
gB = ggplotGrob(carbon_intensity)
gC = ggplotGrob(energy_gdp)
gD = ggplotGrob(carbon_intensity_gdp)
png("relative_dematerial.png", width = 7.5, height = 3.2,  units = 'in', res = 600)
grid.arrange( rbind(
cbind(gA, gB, size = "last"),
cbind(gC, gD, size = "last"),
size = "last")
)
dev.off()
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
# regression
x = log(d$NV.SRV.TETC.ZS)
y = log(d$EG.EGY.PRIM.PP.KD)
r = lm(y ~ x)
r.square = summary(r)$r.squared
x = seq(log(1), log(110), length.out = 50)
x = data.frame(x)
regression = exp(predict(r, x,  interval = "prediction", level = 0.99))
regression = data.frame(x = exp(x), regression)
intensity.r2 =  paste("R^2 == ", round(r.square,2))
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial_value_added.R')
library(data.table)
library(ggplot2)
library(here)
library(gridExtra)
library(grid)
library(gtable)
dir = here()
wd = gsub("figures", "data/world bank", dir)
setwd(wd)
d = fread("data_format.csv")
d = d[country != "Lesotho"]
d$EG.USE.PCAP.KG.OE = d$EG.USE.PCAP.KG.OE*4.187/100
text.size = 10
source('~/Desktop/service/figures/absolute_dematerial_value_added.R')
source('~/Desktop/service/figures/non_production.R')
source('~/Desktop/service/figures/non_production.R')
library(data.table)
library(ggplot2)
library(here)
library(gridExtra)
library(grid)
library(gtable)
library(reshape2)
text.size = 10
# EMPLOYMENT
# energy intensity vs. service employment
###################################################
dir = here()
wd = gsub("figures", "data/world bank", dir)
setwd(wd)
d = fread("data_format.csv")
d$EG.USE.COMM.FO.ZS
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
x = log(d$SL.SRV.EMPL.ZS)
y = log(d$EG.EGY.PRIM.PP.KD*d$EG.USE.COMM.FO.ZS/100)
r = lm(y ~ x)
plot(x, y, log = "xy")
x
y
summary(x)
summary(y)
summary(d$EG.USE.COMM.FO.ZS)
x = x[y != -Inf]
x = x[y != -Inf]
y = y[y != -Inf]
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial_value_added.R')
source('~/Desktop/service/figures/relative_dematerial_value_added.R')
intensity.r2 =  paste("R^2 == ", round(r.square,3))
library(data.table)
library(ggplot2)
library(here)
library(gridExtra)
library(grid)
library(gtable)
library(reshape2)
text.size = 10
dir = here()
wd = gsub("figures", "data/world bank", dir)
setwd(wd)
d = fread("data_format.csv")
# VALUE ADDED
# energy intensity vs. service gdp
###################################################
# regression
x = log(d$NV.SRV.TETC.ZS)
y = log(d$EG.EGY.PRIM.PP.KD*d$EG.USE.COMM.FO.ZS/100)
x = x[y != -Inf]
y = y[y != -Inf]
r = lm(y ~ x)
r.square = summary(r)$r.squared
x = seq(log(1), log(110), length.out = 50)
x = data.frame(x)
regression = exp(predict(r, x,  interval = "prediction", level = 0.99))
regression = data.frame(x = exp(x), regression)
intensity.r2 =  paste("R^2 == ", round(r.square,2))
summary(regression)
summary(r)
source('~/Desktop/service/figures/relative_dematerial_value_added.R')
source('~/Desktop/service/figures/relative_dematerial_value_added.R')
source('~/Desktop/service/figures/relative_dematerial_value_added.R')
source('~/Desktop/service/figures/relative_dematerial_value_added.R')
source('~/Desktop/service/figures/absolute_dematerial_value_added.R')
source('~/Desktop/service/figures/absolute_dematerial_value_added.R')
source('~/Desktop/service/figures/absolute_dematerial_value_added.R')
source('~/Desktop/service/figures/absolute_dematerial_value_added.R')
source('~/Desktop/service/figures/absolute_dematerial_value_added.R')
source('~/Desktop/service/figures/absolute_dematerial_value_added.R')
source('~/Desktop/service/figures/absolute_dematerial_value_added.R')
source('~/Desktop/service/figures/absolute_dematerial_value_added.R')
library(data.table)
library(ggplot2)
library(here)
library(gridExtra)
library(grid)
library(gtable)
dir = here()
wd = gsub("figures", "data/world bank", dir)
setwd(wd)
d = fread("data_format.csv")
text.size = 10
d = fread("data_format.csv")
d$EG.USE.PCAP.KG.OE = d$EG.USE.PCAP.KG.OE*4.187/100
x = log(d$SL.SRV.EMPL.ZS)
y = log(d$EG.USE.PCAP.KG.OE)
r = lm(y ~ x)
r.square = summary(r)$r.squared
x = seq(log(1), log(100), length.out = 50)
x = data.frame(x)
regression = exp(predict(r, x,  interval = "prediction", level = 0.99))
regression = data.frame(x = exp(x), regression)
library(data.table)
library(ggplot2)
library(here)
library(gridExtra)
library(grid)
library(gtable)
dir = here()
wd = gsub("figures", "data/world bank", dir)
setwd(wd)
d = fread("data_format.csv")
d$EG.USE.PCAP.KG.OE = d$EG.USE.PCAP.KG.OE*4.187/100
text.size = 10
# fossil fuel fraction
#############################################################
# regression
x = log(d$SL.SRV.EMPL.ZS)
y = log(d$EG.USE.COMM.FO.ZS)
x = x[d$EG.USE.COMM.FO.ZS > 0]
y = y[d$EG.USE.COMM.FO.ZS > 0]
r = lm(y ~ x)
r.square = summary(r)$r.squared
x = seq(log(1), log(100), length.out = 50)
x = data.frame(x)
regression = exp(predict(r, x,  interval = "prediction", level = 0.99))
regression = data.frame(x = exp(x), regression)
fossil.r2 =  paste("R^2 == ", round(r.square,2))
fossil = ggplot() +
geom_ribbon(data = regression, aes( x = x, ymin = lwr, ymax = upr), fill = "black", alpha = 0.05) +
geom_line(data = regression, aes(x = x, y = fit), col = "grey50") +
geom_point(data = d, aes(x = SL.SRV.EMPL.ZS, y = EG.USE.COMM.FO.ZS), size = 0.1, col = "dodgerblue4") +
scale_x_log10("Service Sector (% of Employment)", breaks = c(2, 5, 10, 20, 50, 100)) +
scale_y_log10("Fossil Fuels (% Total Energy)",
breaks = c(2, 5, 10, 20, 50, 100) ) +
coord_cartesian(xlim = c(4, 100), ylim = c(1.5, 100)) +
ggtitle("D.  Energy from Fossil Fuels") +
theme_bw() +
theme(panel.border = element_rect(color = "black"),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
legend.position = "none",
plot.title = element_text(face="bold", size = rel(1), hjust = 0.5),
axis.line = element_line(color = "black"),
axis.title.x=element_text(vjust=-0.3, size=rel(0.9)),
axis.title.y=element_text(vjust= 1.1, size=rel(0.9)),
axis.text.x = element_text(margin=margin(5,5,0,0,"pt")),
axis.text.y = element_text(margin=margin(3,5,0,3,"pt")),
axis.ticks.length = unit(-0.7, "mm"),
text=element_text(size = text.size, family="Times")) +
annotate("text", x = 5, y = 100, label = fossil.r2, parse = T,family="Times", size=3)
# total energy
####################################################################
# regression
x = log(d$SL.SRV.EMPL.ZS)
y = log(d$EG.USE.PCAP.KG.OE)
r = lm(y ~ x)
r.square = summary(r)$r.squared
x = seq(log(1), log(100), length.out = 50)
x = data.frame(x)
regression = exp(predict(r, x,  interval = "prediction", level = 0.99))
regression = data.frame(x = exp(x), regression)
energy.r2 =  paste("R^2 == ", round(r.square,2))
energy = ggplot() +
geom_ribbon(data = regression, aes( x = x, ymin = lwr, ymax = upr), fill = "black", alpha = 0.05) +
geom_line(data = regression, aes(x = x, y = fit), col = "grey50") +
geom_point(data = d, aes(x = EG.USE.PCAP.KG.OE, y = SL.SRV.EMPL.ZS), size = 0.1, col = "dodgerblue4") +
scale_x_log10("Energy Use per Capita (GJ)", breaks =  c(1,10, 100, 1000)  ) +
scale_y_log10("Agriculture (% Employment)", breaks = c(1, 2, 5, 10, 20, 50, 100)) +
coord_cartesian(xlim = c(1, 1000), ylim = c(5, 100)) +
ggtitle("C.  Services") +
theme_bw() +
theme(panel.border = element_rect(color = "black"),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
legend.position = "none",
plot.title = element_text(face="bold", size = rel(1), hjust = 0.5),
axis.line = element_line(color = "black"),
axis.title.x=element_text(vjust=-0.3, size=rel(0.9)),
axis.title.y=element_text(vjust= 1.1, size=rel(0.9)),
axis.text.x = element_text(margin=margin(5,5,0,0,"pt")),
axis.text.y = element_text(margin=margin(3,5,0,3,"pt")),
axis.ticks.length = unit(-0.7, "mm"),
text=element_text(size = text.size, family="Times")) +
annotate("text", x = 5, y = 100, label = energy.r2, parse = T,family="Times", size=3)
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
blank = ggplotGrob(blank)
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/disag.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/relative_dematerial.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
blank = ggplot() +
geom_blank() +
theme_bw() +
theme(panel.border = element_blank())
source('~/Desktop/service/figures/absolute_dematerial.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
source('~/Desktop/service/figures/absolute_dematerial.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/us_intensity_plot.R')
source('~/Desktop/service/figures/relative_dematerial.R')
library(data.table)
library(ggplot2)
library(here)
library(gridExtra)
library(grid)
library(gtable)
library(reshape2)
text.size = 10
# US sectors by value
#######################################################################
dir = here()
wd = gsub("figures", "/data/bea", dir)
setwd(wd)
d = fread("sector_energy_intensity_va.csv")
boxplot(d$industry/d$service)
d = fread("sector_energy_intensity_pep.csv")
long = data.table(intensity = c(d$industry, d$service),
sector = rep(c("US Industry", "US Services"), each = length(d$year))
)
boxplot(d$industry/d$service)
